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1.0 SUMMARY 


A general procedure is developed to determine the attitude of the longitudinal axis (spin axis) 
and a lateral axis (experiment axis) of a rotating vehicle from the output of solar sensors and a 
lateral magnetometer. 

The orientation of the solar vector with respect to a rocket-axis system is obtained from solar 
sensors located both on the side and on the nose of the rocket. The lateral magnetometer is used 
to obtain the location of the component of the magnetic field vector along the lateral plane of the 
rocket. The solar vector is located in a horizon coordinate system from solar data obtained from 
"The American Ephemeris," and rocket radar inputs. The geomagnetic field vector in the horizon 
coordinate system is calculated from the Jensen and Cain field model using rocket radar data. The 
location of the longitudinal axis and a lateral (experiment) axis of the vehicle in both horizon co- 
ordinates and space fixed coordinates is calculated by the use of various coordinate transformations 
as discussed in references (3) and (4). 


The method is programmed on an IBM 7094 computer. An explanation of the computer opera- 
tions and a copy of the computer program is included in the report. A particular rocket is used as 
a sample case and the results are presented as time histories of the zenith and azimuth angles 
(horizon coordinates) and the Right Ascension and Declination (space fixed coordinates) of both the 
longitudinal and experiment axes. VsiSs / 

^ / 


2.0 INTRODUCTION 



In order to translate and make use of the scientific information obtained from a rocket flight, 
it is often important that the attitude of the vehicle be known at a given time in flight. Various 
instruments can be included in the telemetry package to obtain this information such as; solar 
sensors, magnetometers, horizon sensors, and gyroscopes. This method utilizes solar sensors 
and a lateral magnetometer. 


The solar sensor determines two angles which define two perpendicular planes each containing 
the sun. The intersection of the two planes defines a line pointing at the sun. The field of view of 
the sensor has the shape of an inverted pyramid with a typical angle between the opposite faces 
being 128 degrees. A total of 5 such sensors are necessary for a 4 Pi steradian coverage (see 
Figure 9). The lateral magnetometer simply measures the voltage output which is in turn propor- 
tional to the component of the magnetic field along that particular lateral axis. 

It is necessary to locate at least two spatial vectors with respect to both the rocket and earth/ 
space in order to obtain the orientation of the rocket with respect to earth/space. One vector is 
not sufficient as the results obtained are the same when the rocket is rotated about that vector. 


The two vectors chosen in this procedure are the solar vector and the magnetic field vector. 
The method is therefore applicable only when the sun is within view of at least one solar sensor 
and the lateral component of the magnetic field is monitored by a lateral magnetometer. 


3.0 COORDINATE SYSTEMS 


It is necessary to describe the various coordinate systems referred to in the report. The three 
orthogonal unit vectors in each of the following right-handed systems are defined along the appro- 
priate axes. 

3.1 Rocket-axis Coordinate System (Figure 1) 

T - experiment axis (an axis in the lateral rocket plane chosen as the reference axis 

y - y axis (an axis in the lateral plane orthogonal to the experiment axis and the rocket longi- 
tudinal axis) 

k - rocket longitudinal axis (spin axis of vehicle) 

3.2 Horizon Coordinate System (Figures 5) 

T - East axis (an axis in the local horizon plane in the east direction) 

J - North axis (an axis in the local horizon plane in the north direction) 

k - Vertical axis (an axis perpendicular to the horizon plane which is positive when measured 
upward) 

3.3 Equatorial Coordinate System (Figure 5) 

T - Equator axis (an axis along the equatorial plane perpendicular to the North Celestial Pole 
axis at the local meridian of the rocket.) 

J - East axis (an axis in the equatorial plane orthogonal to the North Celestial Pole axis and 
the equator axis). 

k - North Celestial Pole axis (an axis perpendicular to the equatorial plane measured from the 
center of the earth to the North Celestial Pole). 

3.4 Space- Fixed Coordinate System (Figure 8) 

j - Vernal equinox (an axis measured from the center of the earth towards the vernal equinox). 

y — y * (an axis on the equitorial plane orthogonal to both the Vernal equinox and the North 
Celestial Pole axes). 

k - North Celestial Pole axis (same as in 3.3). 

4.0 SYMBOLS 

A Azimuth angle 

Aj Lateral solar sensor angle (Case I sensor) 

A 2 Longitudinal solar sensor angle (Case II sensor) 
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B t Lateral solar sensor angle (Case I sensor) 

B 2 Longitudinal solar sensor angle (Case II sensor) 

E Experiment axis vector 
H Magnetic field vector 

T, I, k Units vectors in right handed orthogonal coordinate systems 
La Rocket latitude 
Lo Rocket longitude 
r Rocket longitudinal axis vector 
s Solar vector 
t Greenwich Sidereal time 

c 

V Arbitrary vector 

Xj Component of the solar vector along the experiment axis 

X 2 Component of the magnetic field vector along the experiment axis 

X 3 Component of the solar vector along the east axis 

X 4 Component of the magnetic field vector along the east axis 

x' Coordinate axis in lateral rocket plane 

Yj Component of the solar vector along the y axis 

Y 2 Component of the magnetic field vector along the y axis 

y 3 Component of the solar vector along the north axis 

Y 4 Component of the magnetic field vector along the north axis 

y Coordinate axis orthogonal to the Experiment axis and the rocket longitudinal axis 
y' Coordinate axis orthogonal to the x' axis and the rocket longitudinal axis 
z Zenith angle 

Z, Component of the solar vector along the rocket longitudinal axis 
Z 2 Component of the magnetic field vector along the rocket longitudinal axis 
Z 3 Component of the solar vector along the North Celestial Pole axis 
Z 4 Component of the magnetic field vector along the North Celestial Pole axis 
. Right Ascension 
Declination 
Hour angle 

Azimuthal angle in the rocket lateral plane measured from the experiment axis 

/' Azimuthal angle in the rocket lateral plane measured from the x' axis (Case II 
sensor only) 
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<p si Angle between experiment and x' axis (Case I sensor) 

4> S2 Angle between experiment and x' axis (Case n sensor) 
o Elevation angle from rocket lateral plane 
All angles are in radians unless otherwise specified. 

Subscripts 

E Pertaining to the experiment axis vector 
M Pertaining to the magnetic field vector 
R Pertaining to the rocket longitudinal axis vector 
S Pertaining to the solar vector 
V Pertaining to an arbitrary vector 

ER Pertaining to both the experiment axis vector and the rocket longitudinal axis vector 
ZA Pertaining to a vector in the horizon coordinate system 
eb Pertaining to a vector in the equatorial coordinate system 

5.0 METHOD OF ANALYSES 

5.1 Solar Vector-Horizon Coordinate System (Figure 5) 

The position of the solar vector with respect to the horizon coordinate system can be found 
from solar data in reference (1) and the rocket radar inputs; latitude, longitude, and altitude as a 
function of time. 

5.2 Magnetic Field Vector — Horizon Coordinate System (Figure 6) 

The magnetic field vector can be defined in the horizon coordinate system by use of the spherical 
harmonic analysis in reference (2). For a given latitude, longitude and altitude, the three earth- 
fixed components of the magnetic field vector can be obtained. 

5.3 Solar Vector - Rocket-Axis System (Figures 1 and 2) 

Each solar sensor measures the angle of incident sunlight with respect to the rocket-axis co- 
ordinate system along two perpendicular planes. For the sensor mounted along the side of the 
rocket, these planes will be the lateral plane and a plane containing the sensor and the vehicle longi- 
tudinal axes (Figure 1). For the sensor mounted on the nose of the vehicle, each plane contains the 
longitudinal axis and an arbitrary axis in the lateral plane. The lateral axes for the two planes 
are mutually perpendicular (Figure 2). Thus the solar vector can be totally defined with respect 
to the rocket-axis system by either type of sensor as long as the vector is within the viewing range 
of the sensor. 
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5.4 Magnetic Field Vector - Rocket-Axis System (Figure 3) 


The lateral magnetometer measures the component of the external magnetic field along the lateral 
axis. The output of the magnetometer is sinusoidal for a rolling vehicle with the maximum or mini- 
mum corresponding to the alignment of the magnetometer along or opposite the projection of the 
magnetic field vector on the lateral plane. At this point in the analysis, only the angular location 
of the magnetic field vector in the rocket lateral plane (</> ) is known. The elevation of this vector 
( f u ) must be determined to totally define it in the rocket axis system. The magnetic field vector 
can then be totally defined in the rocket- axis system by equating the dot product of the solar and 
magnetic field vectors in the rocket-axis system with the known value of the dot product of the two 
vectors in the horizon coordinate system and subsequently obtaining the elevation angle of the mag- 
netic field vector ( ). 

5.5 Longitudinal and Experiment Axes — Horizon Coordinate System (Figure 5) 

Two spatial vectors are at this point completely defined in the rocket-axis and horizon coordi- 
nate system. It is then necessary to obtain each rocket-axis unit vector in terms of horizon unit 
vectors. To do this, a new coordinate system is defined utilizing functions of the solar and mag- 
netic field vectors as unit vectors. Each rocket coordinate axis can then be represented vectorially 
in this new coordinate system and each unit vector of this system is known in terms of the horizon 
system unit vectors. Thus each rocket-axis coordinate can be defined in the horizon coordinate 
system. 

5.6 Longitudinal and Experiment Axes — Space Fixed Coordinates (Figures 5 and 8) 

Each rocket system coordinate axis can be determined in space fixed coordinates by first con- 
verting to the equatorial coordinate system. This is obtained by utilizing the relationship between 
the unit vectors of the horizon and equatorial coordinate systems. A simple conversion can then be 
made to obtain the results in space fixed coordinates (right ascension and declination) from equa- 
torial coordinates (hour angle and declination). 

6.0 ANALYSIS 

6.1 Solar Vector — Rocket-Axis System 

The equations for obtaining the position of the solar vector in the rocket-axis coordinate system 
can be divided into the two following possible cases: 

Case I. The solar sensor senses along an axis normal to the vehicle longitudinal axes (figure 1 ). 
Case II. The sensor senses along the longitudinal axis (figure 2). 

Case I. / s = i sl i B, (1) 

The solar vector in a coordinate system consisting of the sensor axis (i ), rocket longitudinal 
axis (k ) and a third orthogonal axis (I ) is: 
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(2) 


S = cos £) s cos Bj i + cos sin Bj j + sin d g k 
The component of the vector in the rocket-sensor plane is then: 

_ cos d g cos B 1 i + s in d g k 

k — > 

y cos 2 d g cos 2 Bj + s in 2 d g 

Equation (3) is then dotted with the unit vector T to obtain cos A t 

__ _ cos d g cos Bj 

cos A, = S ik • T = r 

ycos 2 ^g cos 2 Bj + sin 2 d g 

therefore 


tan 2 d 

see 2 Aj = 1 + 

cos 2 Bj 


and 


d g tan -1 (cos tan Aj) 


(3) 


(4) 


( 5 ) 

( 6 ) 


Case II. The solar vector in the coordinate system consisting of the rocket axis (k ), which is also 
the sensor axis, and the two orthogonal axes (T, J) along which the incident sunlight is measured 
is: 


S - cos fV s cos & i + cos sin <£' s j + sin d g k 

The component of the vector in the rocket, x' plane is: 


cos d s cos 4 ' g i + sindgk 
]/~cos 2 d cos 2 + sin 2 d s 


Equation 8 is then dotted with the unit vector k 


_ _ sin d 

cos (-A 2 ) : S ik - k = - — 

y cos 2 d s cos 2 <^g + sin 2 d s 


(7) 


( 8 ) 


(9) 
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then 


tan 2 (- cot 2 6^ cos 2 4>' s 

The component in the rocket, y' plane is: 


cos 0^ sin <£g j + sin k 
(cos 2 # s sin 2 4>' s + sin 2 ^ s ) l/2 


( 10 ) 


( 11 ) 


Equation (11) is then dotted with k 


sin fL 


cos B 2 - * k 


(cos 2 6L sin 2 4' i sin 2 0 ) l y 2 


( 12 ) 


then 


tan 2 B 2 cot 2 (J sin 2 (/;* 


(13) 


,f can be solved by combining equations (10) and (13) 


-tan B 

'/■' tan' 1 ' 2 

‘tan 


(14) 


where: 


i • 

1 s 



if: A 2 

< 0, 

B 2 

> 0 


180° 4 

/ • 

7 s 

A 2 

> 0, 

B 2 

> 0 


180° i 

'/*8 

A 2 

> 0, 

B 2 

< 0 

‘/'s ^ 

360° 4 

1 S 

A 2 

< o, 

B 2 

< 0 


then 


S ^ S ^S2 


(15) 


Equation (14) is substituted into (13) to obtain: 


°s * an ' 


t an 2 A 2 + t an 2 B 2 


(16) 
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The solar vector in the x , y, z coordinate system (figures 1 and 2) is: 


® ~ ^1 ^EXPERIMENT + ^1 + ^SPIN 

where: 

Xj = cos 0 S cos t/> s 
Y. - cos 0 sin (p 

1 S *S 

Z x □ sin # s 

o s and 0 s are obtained from equations (1) and (6) or (15) and (16) depending on whether case I or 
case II is applicable. 

6.2 Magnetic Field Vector - Rocket- Axis System (Figure 3) 

When the magnetic vector is aligned along the maximum projection of the magnetic vector on 
the lateral plane, a maximum voltage output is telemetered to the ground station. At this time, the 
magnetometer is sensing at which is the angle between the magnetometer axis and the refer - 

max 

ence axis. Conversely, when it is aligned opposite to this projection, a minimum output will result 
corresponding to .A time history of ch can then be obtained by interpolation between the 

i n M 

maximum and minimum values. 

The magnetic vector in the rocket axis system is: 

H “ ^2 ^EXPERIMENT 4 ^2 ^ Y + ^2 ^SPIN (18) 

where: 

X 2 = cos o u cos 0^ 

Y 2 cos U M sin c/^ 

Z 2 = sin 0 U 

with o yet to be determined. 

6.3 Solar Vector - Horizon Coordinates (Figures 4 and 5) 

The hour angle of the solar vector is obtained by noting figure 4. 

( S = a S - f G - L ° (19) 
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where <x s and t G are found from reference (1) using the appropriate universal flight time and Lo 
is obtained from the rocket radar data. 

Figure 5 is a representation of any vector, v , in both a horizon and an equatorial coordinate 
system. It can be written as: 

V ZA = sin Z sin A I East + sin Z cos A T North + cos Z k Verticfll (horizon) (20a) 

V 8 = cos £ cos S i Equator + sine cos 8 I East + sin 8 k North Polc (equatorial) (20b) 


Each unit vector in one coordinate system can be represented by the unit vector in the other: 


^East ^East 


j M - — sin La i_ A + cos La k__ _ _ , 

J North Equator North Pole 


(21a) 


k„ _ . - cos La i-. ^ + sin La k M . „ , 

Vertical Equator North Pole 


i La i. 


eoc T o 


^ Eas t 


(21b) 


k_, . _ , - cos La i_, + sin La k., ^ . . 

North Pole J North Vertical 


If the vector in question is the solar vector, the horizon components can be found by taking the dot 
product of equation (20b) and the unit vectors in (21a). 

Then 


S - X 3 i £ast + ^3 ^fj or th * ^3 ^Vertical 


where: 

X 3 " g ^‘ r East " Sin f S C ° S S S 

Y 3 ~ ^rI North = - Sin La COS 6 S COS S + COS La Sin 

Z, = S « *k„ , cos La cos £ c cos b- 4 sin La sin 

3 Vertical =» S a 

and La is obtained from the rocket radar data. 


( 22 ) 
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6.4 Magnetic Field Vector - Horizon Coordinates (Figure 6) 

The components of the magnetic field vector in the horizon coordinate system (x 4 , Y 4 , z 4 ) are 
calculated using the Jensen and Cain, 48 term Expansion in Spherical Harmonics of the Magnetic 
Field in 1960 (reference 2). These equations have been programmed on the 7094 digital computer 
and are included as a subprogram of the general attitude program. The inputs to this subprogram 
are the vehicle's latitude, longitude, and altitude which are obtained from the rocket radar data. 

The Magnetic field vector is then: 


^ ^4 *East ^4 ^North ^4 ^Vertical 


(23) 


6.5 Elevation of Magnetic Field Vector - Rocket- Axis System (Figure 3) 

The dot product of two vectors is the same in any coordinate system; that is, 

/C . _ /Q . tl\ 

' ' Ro c k e t Coordinates ~ V 'Horizon Coordinates 

Thus dotting equation (17) with (18), and(22) with (23) the following is obtained: 


X 1 X 2+ Y 1 Y 2+ Z 1 Z 2= X 3 X 4+ Y 3 Y 4+ Z 3 Z 4 


or 


cos 0 U cos <£ m Xj + cos e H sin ^ Yj + sin 0^ Z x = X 3 X 4 + Y 3 Y 4 + Z 3 Z 4 


(24) 


where : 


^ is the only remaining unknown. 
Equation (24) is solved for sin e u and cos 



n 4A m C im) 1/2 

S1 M “ 2^ 

(25a) 

where: 

n _ ” B 2M ± ^ B 2M “ 4A M C 2M^ 1/2 

C ° M_ 2A, 

(25b) 


\ = Z? + ( X ! cos <4, + Yj sin <^) 2 



B 1(J = - 2(X 3 X 4 + Y 3 Y 4 + Z 3 Z 4 ) Z, 
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B 2M = - 2 < X 3 X 4 + Y 3 Y 4 + Z 3 Z 4> <*i COS <*W + Y , sin 

C 1U = ( X 3 X 4 < Y 3 Y 4 + Z 3 Z 4 ) 2 " ( X ! COS ^ 4 Yj sin ^> 2 

C 2M =< X 3 X 4^ Y 3 Y 4^ W“ Z 1 2 

Solving equation (25b) results in two solutions for cos o u . Since, by definition, 0 U must be 

either in the first or fourth quadrant (figure 3), cos 6 U must be positive. It is possible for both 

values of cos o u to be positive corresponding to a uniformly rotating case and an extreme coning 
case. The value of cos a u corresponding to the extreme coning case can easily be eliminated by 
observing the continuity of the results. Once the correct value of cos o u is known, the correspond- 
ing value of sin t? u from (25a) can be obtained. 

6.6 Rocket Longitudinal and Experiment Vectors - Horizon Coordinates (Figures 5 and 7) 

In order to obtain the Zenith and Azimuth angles of the rocket and experiment vectors, a new 
coordinate system is defined consisting of the unit vectors 


S X H 


Hx(SxH) 

(1 _ (S-H) 2 ) 1/2 _ 


(1 _ (S -H ) 2 ) 1 /2 


These unit vectors can be represented in both the rocket-axis system and the horizon coordi- 
nate system. Thus the direction cosines of any arbitrary vector will be the same in either coordi- 
nate system. From figure 7, any vector in this new coordinate system becomes: 


V (V *H) [H] 


V • (S' x H) 

S x H 

(V-S) -(S-H) (V * H) 

i , 

S - (S‘H)H 

(1 - (S -H ) 2 > 1/2 

(1 - (S-H ) 2 ) 1/2 

(1 - (S-H ) 2 ) 1/2 

(1 -(S-H ) 2 ) 1/2 


( 26 ) 


The cross product of equations (17) and (18) yields: 


S x H 


<V a - W 


Expe r imrn t 


(ZjXj - XjZ 2 ) j Y + (X,Y 2 


•w 


(27) 


The rocket longitudinal and experiment vectors are then dotted with equations (17), (18), and (27) 


RS Zj 

R H r Z 2 ( 28 a) 

R-(S . H) - X,Y 2 - Y,X 2 

ES Xj 

EH = X 2 ( 28 b) 

E-(S > H) =Y,Z 2 -Z 2 Y 2 
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Equations (28a) and (28b) are substituted into equation (26) to obtain the rocket longitudinal and 
experiment vectors. 


X.Y, - Y X 

R - 7 [HI , 12 12 

S X H 

Z, - (S-H) Z 2 

S - (S-H) H 

(1 - (S-H) 2 )'/ 2 

(1 - (S-H) 2 ) 1/2 

(1 - (S-H) 2 )*/ 2 

(1 - (S-H) 2 ) 1/2 


Y, Z - Z Y 

E - X [H] + 

Sx H 

X, - (S-H) X 2 

S - (S-H) H 

(1 - (S-H) 2 ) 1/2 

(1 _ (S-H) 2 ) 1/2 

(1 - (S-H) 2 ) 1/2 

(1 - (S-H) 2 ) 1/2 


(29a) 


(29b) 


The cross product of equations (22) and (23) yields: 


s*H = (Y 3 z 4 -z 3 Y 4 ) i East 


+ - W Wu, 


+ (X 3 Y 4 


Y,xj k 

3 4 Verti 


(30) 


The components of the rocket longitudinal and experiment vectors in the horizon coordinate system 
can be obtained by dotting each vector (as denoted by both equations (29a), (29b) and equation (20a)) 
into each horizon unit vector using equations (22), (23), and (30) 


R-i 


East 


sin Z„ sin A„ - Z X, + 
R R 2 ^ 


( X . Y 3 -W ( y 3 Z 4 - Z 3 V 
1 - (S-H) 2 


+ 


[Z, - (S-H) Z 2 1 (X 3 - (S-H) X 4 ] 
1 - (S-H) 2 


(31a) 


R --> Nort h = s in Z R cos A r - Z 2 Y 4 + 


(XiY 2 -Y,X 2 ) (Z 3 X4-X 3 Z 4 ) 
1 - (S-H) 2 


[2T 1 _ (S-H) Z 2 ] LY 3 - (S-H) Y 4 ] 


1 - (S-H) 2 


(31b) 


IMc . 1 ; COS Z Z + 

Vertical R 2 4 


< X , Y 2- W (X 3 Y , "W 

1 - (S-H) 2 


fZj - (S-H) Z 2 J [Z 3 - (S-H) Z 4 
’ 1 - (S-H) 2 


(31c) 
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! 

* 7 . „ . , v v (V i Z 2 -Z 1 Y 2 )(Y 3 7 4 -Z 3 Y 4 ) 

E • *£»*, T Sln Z E s,n *£ = X 2 X 4 + r 

1 - (S • H) 2 


i 

! . 

1 X, - (S • H) Xj] [X 3 - (S • H) X 4 ] 
1 - (S • H) 2 

(32a) 

* * 

5 - (Y.Z 2 -Z.Y,)(Z,X 4 -X,Z 4 ) 

® * 1 North " S * n ^E COS Ag — X2 Y 4 + 

1 - (S • H) 2 


! 

[x, - (S - H) X 2 ] [Y 3 -(S-H)Y 4 ] 
1 - (S • H) 2 

(32b) 

1 

! 

| 

__„ 7 V7 ( Y i Z 2- Z .V 2 )(X3Y 4 -Y 3 X 4 ) 

E k v. „ical ' COS Z E - X 2 Z 4 + r-3 

1 - (S • H) 2 


| 

[X,-(S-H)X 2 ] [Z 3 - (S • H) Z 4 ] 
1 - (S • H) 2 

where : 

(32c) 

! 

S - H X 3 X 4 + Y 3 Y 4 + z 3 z 4 


1 

From equations (31c) and (32c), the Zenith angle of the experiment or rocket longitudinal axis is: 

i 

; 

i 

z., c - tan = 

L cos Z ER 

! 

(33) 

1 

- 

From Figure 5, if: 



^ER > 0 ^ER = ^ER 



^ER ® ^ER 180 4 Z RR 


! 

Equation (33) is substituted into (31b) and (31c) for the rocket longitudinal vector 
(32c) for the experiment vector to obtain: 

or (32b) and 
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( 34 ) 


cos A er 


(sin 


Z £R cos A £R ) 
sin Z ER 


(35) 


The Azimuth of the rocket or experiment vector is obtained by dividing equation (34) by (35) 


A 


ER 


tan 


sin A £R 
cos A £R 


= tan 1 (tan A ER ) 


(36) 


tan A £R > 0, sin A £R > 0 
tan A er < 0, sin A £R > 0 
tan A £R > 0, sin A £R < 0 
tan A £R < 0, sin A £R < 0 


then: A £R - A £R 

^ER “ + ^ER 

A £r =180 + A £R 
^er = + \r 


6.7 Rocket Longitudinal and Experiment Vectors - Equatorial Coordinates (Figures 4 and 5) 

By taki ng the dot product of the rocket longitudinal and experiment vectors in horizon coordi- 
nates (equation 20a) and each equatorial unit vector (equation 21b), the components of the experiment 
and rocket longitudinal vectors in equatorial coordinates can be obtained. 


V ER * ‘‘‘Equator 


cos k. 


- sin La sin Z FP cos A F 


+ cos La cos Z ER 


^ER * J Eas t 


V • w 

V ER “-North Pole 


= sin e £R cos S £R = sin Z £R sin A £R 


- sin 8 £R - cos La sin Z ER COS A ER 
+ sin La cos Z £R 


(37a) 

(37b) 

(37c) 
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The declination of the experiment and rocket longitudinal vector can be obtained from equation (37c) 


S ER " tan * 


sin ii- 


(t - sin 2 Sg ,) 172 


where -90 c ' < i> ER < 90° 


(38) 


The hour angle can be found by substituting equation (38) into (37a) and (37b). 


sin e. 


(cos sin « ER ) 


(39) 


E R 


(cos S ER cos € er ) 
cos 


(40) 


Therefore dividing equation (39) by (40) 


sin t. 


= tan J (tan t ) 


ER 


if: 


ER 


“ ER 


ER 


ER 


> 0 

then: < ER 

* ER 


> 0 

ER ~ 

180 ° 

+ € 

< 0 

l er 

180° 

+ f 

< 0 

f ER " 

360 " 

+ £, 


(41) 


The Right Ascension can then be obtained from figure 4. 


*er ' € ER + t G t ho 


(42) 


7.0 COMPUTER PROGRAMMING 

The procedure has been written in a program for an IBM 7094 computer in Fortran IV. It 
could also be applied to other computers with only a few format changes. A copy of the computer 
program is shown below along with an explanation of the data inputed in the particular sample case, 
the results of which are shown in Appendix B. 
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Tract Program Inputs 

The Tract Program (Fortran IV), which determines the body orientation of the rocket, has 
been designed to run on the Moonlight System (IBM 7040-7094 direct Couple system). 

Two BCD tapes are inputed to the program and are set up on tape units 18 and 19. 

I. Input Tapes 

A. Tape 18 contains the telemetry data and each record should contain time, latitude, longi- 
tude, and altitude in this order. The units are secs, (after launch), degrees, and feet, re- 
spectively. 

B. Tape 19 contains the sunsensor data and each record should contain time (after launch), 
the eye or sensor number which was viewing at this time (1, 2, 3, or 4), the lateral solar 
sensor angle, and the longitudinal solar sensor angle (see Figures 1 and 2). The units are 
secs, none (integer), degrees, and degrees, respectively. 

II. Input Cards 

A. CARD 1 - Format (215, 7F10.5) Nine values are to be read in. 

(1) if this value is nonzero the program will go through a Calcomp plotting sequence 
yielding an output tape at 200 B.P.I. which can be plotted on an uncoupled Calcomp 
Plotter, 

(2) this value determines the number of cards to be read in section B, 

(3) the initial time where the program will start computing, 

(4) the final time the program will execute, 

(5) the constant time increment for computation (the three above values are in secs), 

(6) the angle between the experiment and ! x T axis, Case I sensor, Eye #1, 

(7) the angle between the experiment and ’x f axis, Case I sensor, Eye #2, 

(8) the angle between the experiment and f x' axis, Case I sensor, Eye #3, 

(9) the angle between the experiment and f x T axis, Case II sensor (the four above values 

are in degrees). 

B. The number of cards to be read in this section equals the number in position two of Card 1. 
Each card contains two values (2F10.5), time after launch (secs), and the angle between the 
experiment axis and the maximum or minimum (alternates) projection of the magnetic 
vector (degrees), respectively. 
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C. Seven cards of three numbers each (3F10.5) are read in this section. The times are read 


in hours, minutes, and seconds, respectively, on the card. 
Card 1 - Greenwich sidereal time of launch day. 

Card 2 - Greenwich sidereal time of following day. 

Card 3 - Right ascension of sun for launch day. 

Card 4 - Right ascension of sun for following day. 

Card 5 - Declination of sun for launch day. 

Card 6 - Declination of sun for following day. 

Card 7 - Universal time of launch. 


HI. For Calcomp Plotting the following additional cards should be added. 


Card 1 - 
Card 2 - 


Card 3 - 


Card 4 - 


Card 5 - 


A string of Hollerith characters to be printed as the title at the top of the graph. 
The characters should be centered in the 72 columns of the punched card. 

One integer value specifying the number of plots (15). The following cards 
should be added in sequence for each plot, e.g., if there are three plots there 
should be three sets of the following cards. 

A string of Hollerith characters to be printed as the label for the abscissa of 
the graph. The characters should be centered in the 72 columns of the punched 
card. 

A string of Hollerith characters to be printed as the label for the ordinate of 
the graph. The characters should be centered in the 72 columns of the punched 
card. 

Eight integer values (715); (1) any non zero value, (2) I- Printed values of the 
ordinate are 10.**I times the actual values, (3) J-Printed values of the ordinate 
have J digits following the decimal point, (4) K- Printed values of the abscissa 
are 10.**K times the actual values, (5) M- Printed values of the abscissa have 
M digits following the decimal point, (6) number of horizontal grid lines, (7) 
number of vertical grid lines, (8) NC- number of curves on this plot (maximum 
8 ). 
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Card 6 - Can contain from two to ten integer values (1015) depending upon NC (#8, Card 5); 

(1) integer corresponds to the variable to be plotted as the abscissa for the plot, 

(2) (10) the integer values corresponding to the variable to be plotted as the ordi- 
nate or ordinates for the plot, e.g., if there are four curves on this plot then there 
are four ordinate integer values. 

Note - The variables which can be plotted and their integer designations follow: 

Integer Symbol 

1 
2 

3 

4 

5 

6 

7 

8 
9 

Card 7 - Six values (6F10.5), (1) Spacing between horizontal grid lines, (2) Spacing between 
vertical grid lines (inches), (3) maximum value to be printed for the abscissa, 

(4) minimum values to be printed for the abscissa, (5) maximum value to be printed 
for the ordinate, (6) minimum values to be printed for the ordinate. 

A copy of the computer program follows: 


Variable Computed 
Time (sec) 

Azimuth, Rocket Longitudinal Axis (deg) 
Azimuth, Experiment Axis (deg) 

Zenith, Rocket Longitudinal Axis (deg) 

Zenith, Experiment Axis (deg) 

Right Ascension, Rocket Longitudinal Axis (HRS) 
Right Ascension, Experiment Axis (HRS) 
Declination, Rocket Longitudinal Axis (degrees) 
Declination, Experiment Axis (deg) 
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$ JOB 1091C003 405LAHN 

$SETUP I BFTC 18 TD SAVE REEL OOOOOIBCO 

SSETUP I BFTC 19 TD SAVE REEL OOOOOIBCD 

SEXECUTE I B JOB 

$ I B JOB GO, MAP, SOURCE 

$ I BFTC TRACT NOLI ST ,NOREF , NODECK, M94 

DIMENSION Tl(600 ) , PHI Ml { 600 ) , XHR ( 7 ) , XM IN ( 7 ) ,SEC(7),T (7) ,NSCALE (5), 
XT2 { 100 ), XLAT1 ( 1000 ), XLONG 1 (1000), ALT 1( 1000 ),PHIS1( 1000), 

XTHETAK 1000 ) ,CPL VAR (9 , 1000 ) , AA ( 12 ) , BB ( 12 ) ,CC ( 12 ) , BCDX ( 10 ) , 

XEPS IL ( 2 ) , 

XC0SA2M ( 2 ) , SINAZM ( 2 ) , TANAZM( 2),AZM(2), S I NOEL ( 2 ) »CODS I E ( 2 ) , 

X 10 ( 8 ) , 

XCODCOE ( 2 ) ,DELTA ( 2 ) , COSO EL ( 2 ) , SI NEPS ( 2 ) ,COSEPS ( 2 ) ♦ TANEPS { 2 ) , 

XC0SZEN{2),SIZC0A(2),SIZSIA(2),ZEN(2),SINZEN(2),ALPHA(2),XH(2), 

XXM(2),XS(2),IX(2),NX(2) 

1 FORMAT (F10.2,I10,2F10.2) 

2 FORMAT ( 4E 12. 5 ) 

3 FORMAT (12A6) 

4 FORMAT ( 1015) 

5 FORMAT ( 6F 10 *5 ) 

6 FORMAT ( 2F 10 • 5 ) 

7 FORMAT ( 3F 10 • 5 ) 

8 FORMAT (2I5,7F10#5) 

9 FORMAT ( 10X,4HTIME, 10X, 14HZENITH (ROCKET ) ,10X,15HAZ I MUTH( ROCKET ) , 
X10X, 18HZENITH( EXPERIMENT ) , 10X, 1 9HAZ I MUTH ( EX PER I MENT ) ) 

10 FORMAT (9X,F6. 2, 10X , F7 . 3 , IX , 3HDEG , 14X , F7 .3 , IX , 3HDEG , 15 X , F7 . 3 , IX, 
X3HDEG, 18X,F7.3, 1X,3HDEG) 

15 FORMAT ( 19X»23HRIGHT ASCENSI ON ( ROCKET ) »4X » 19HDECL I NAT I ON (ROCKET ) , 
X4X , 27HR I GHT ASC ENS ION ( EXPER IMENT ) , 4X , 23HDEC L I NAT I ON ( E XPER I MENT ) ) 

16 FORMAT ( 19X,I2,3HHRS,2X,I2,3HMIN,2X,F6.3,3HSEC,8X,F7.3,1X,3HDEG, 

X9X , I 2, 3HHRS, 2X,I2,3HMIN,2X,F6*3» 3HSEC, 12X,F7#3*1X ,3HDEG ) 

17 FORMAT < 10X,8HF0R TI HE , 2X, F6.2, 2X, 61HTHE COSINE OF THE AZMITH IS GR 
XEATER THAN 1* COMPUTATION STOPS) 

19 FORMAT ( 1H 1 ) 

20 FORMAT (1H0) 

21 FORMAT ( 4 1H BOTH ROOTS NO GOOD, POINT SKIPPED, TIME=,F10.5) 

22 FORMAT (41H TWO NEGATIVE ROOTS, POINT SKIPPED, TIME=,F10.5> 

23 FORMAT ( 4 1H TWO POSITIVE ROOTS, POINT SKIPPED, TIME=,F10.5> 

30 FORMAT ( 4 1H RAD21 IS NEGATIVE, POINT SKIPPED, TIME=,F10.5> 

31 FORMAT (4 1H RAD22 IS NEGATIVE, POINT SKIPPED, TIME=,F10.5) 

DATA ( BCDX ( J ) , J=l, 10 ) /0120000000000, 0150000000000 ,0 160000000000 , 

XO 170000000000, 0320000000000, 0350000000000, 0360000000000, 
X0370000000000, 0520000000000 , 0550000000000/ 

NOUNT =0 
NOUNT 1= 1 
NOUNT 2=0 

RE AD (5,8) NPL0T,N1,TIN»TFIN,TI NCRE ,SSA,SSB,SSC,SSD 
RE AD (5,6) (T1(J) ,PHIM1(J),J=1,N1) 

RE AD (5,7 ) (XHR( J) ,XMIN( J),SEC( J) , J=l,7) 

WRITE(6, 19) 

JUMP=0 

C SUNSENSOR ROUTINE 
C TWO CASES ARE POSSIBLE 
C CASE I IS A LATERALLY MOUNTED SENSOR 
C CASE I I IS A LONGITUDINALLY MOUNTED SENSOR 
C0NV=57. 295779 
PI=3. 14159265 
PHI SI A=SSA/CGNV 
PHI S1B=SSB/C0NV 
PHIS1C=SSC/ CONV 
DO 210 J= 1 , 7 

210 T ( J )=XHR( J)*3600.+XMIN< J)*60. + SEC( J) 
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T { 5 ) =T ( 5 ) / 15 • 

T ( 6 ) =T ( 6 ) / 15 • 

11 READ (19,1) TIME, NEY E , C , D 
IF(TIME-TIN)11, 12,12 

13 RE AD ( 19,1) TIME , NEY E , C ? D 

12 TL=T IME 

IF (TL-TFIN ) 14, 14,90 

14 A=C /CONV 
B=D/CONV 

T 1 ( NOUNT 1 ) =TL 
TANA=SIN(A)/COS(A) 

TANB=SIN(B)/CQS(B) 

C TEST FOR CASE NUMBER 

IF (NE YE ♦ EQ .4 ) GO TO 66 
C COMPUTATIONS FOR CASE I 

54 THETAS=ATAN(COS(B)*TANA) 

IF (NEYE-2) 51,52,53 

51 PHI S= B+PH I S 1 A 
GO TO 80 

52 PHI S = PHI S 1 B+B 
GO TO 80 

53 PH I $ = PHI S 1C + B 
GO TO 80 

C COMPUTATIONS FOR CASE II 

66 THET A$=AT AN ( 1 . / SORT ( TANA** 2+TANB**2 ) ) 

PHISP=ATAN(-TANB/TANA) 

IF ( A.LT.O. ) GO TO 67 
PHISP=PI+PHISP 
GO TO 70 

67 IFtB.LT.O. ) GO TO 69 
PHI SP =PH I SP 

GO TO 70 

69 PH ISP=2.*PI+PHI$P 

70 PHIS=PHISP+SSD 

C XlvYlfZl, ARE COMPUTED( SOLAR VECTOR, ROCKET AXIS SYSTEM) 

80 PHI SI ( N0UNT1 ) = PHIS 
THETA 1 (NOUNT 1 )=TH ETAS 
NOUNT 1=N0UNT 1+1 
GO TO 13 

90 NOUNT 1=N0UNT1-1 
N T 2 C N T = 1 
TF IN2=TF IN+. 2 
T IN2 = T IN-* 2 

25 READ (18,2) T 2 1 , X L AT 2 , X L 0NG2 , XA LT 2 
I F ( T I N2-T 21)26,26,25 

26 READ ( 18,2) T2 1 , X L AT2 , X L 0NG2 , XA LT 2 
IF (T21-TFIN2)27,28,28 

27 T2(NT2CNT)=T21 
XLAT1(NT2CNT )=XLAT2 
XL0NG1 (NT2CNT ) =XL0NG2 
ALT1(NT2CNT)=XALT2 

NT 2CNT=NT2CNT + 1 
GO TO 26 

28 N2=NT 2CNT- 1 
98 XNOUNT =N0UNT2 

TD=T IN+XNOUNT*TINCRE 
IF(TD*GT.TFIN) GO TO 800 
NOUNT 2=N0UNT2+1 

CALL TIMC0R(TD,T1, PHI SI, NOUNT 1, PHI S2) 

CALL T IMCOR ( T D, Tl, THET A 1, NOUNT 1, THET A2 ) 

100 X1=C0S(THETA2)*C0S(PHIS2) 

Y1=C0S(THETA2 )*SIN(PHI$2 ) 

2 1 = S I N (THE TA 2 ) 

C MAGNETIC FIELD COMPONENT SECTION 

C THIS SECTION COMPUTES X4,Y4,Z4 WITH LAT I TUDE , L ONG ITUDE , AND ALTITUDE IMPUTS 
C MAGNETIC FIELD VECTOR, EARTH-FIXED AXIS SYSTEM 
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CALL TIMC0R(TD,T2,XLAT1,N2,XLAT) 

CALL TIMCOR(TD,T2»XLONGl,N2, XLONG) 

CALL T IMCOR ( TO* T2» ALT1,N2, ALT ) 

X L ON G=— XLONG 
ALT=ALT/3280.833 

CALL FIELD (XL AT, XLONG, ALT, BT, BP, BR,B2) 

X4=BP/B2 

Y4=-BT/B2 

Z4=BR/B2 

C THIS SECTION LOCATES THE SOLAR VECTOR RELATIVE TO AN EARTH-FIXED 
C AXIS SYSTEM 

C NUMBER =NUMBER OF TIMES THROUGH PROGRAM 

C T ( 2 ) - GREENWICH SI DERI AL TIME OF FOLLOWING DAY (0 HRS. UT ) 

C T(2 ) - GREENWICH SIDERIAL TIME OF FOLLOWING DAY (OHRS. UT ) 

C T ( 3 1 - RIGHT ASCENSION OF SUN FOR LAUNCH DAY (0 HRS. UT J 

C T ( 4 1 - RIGHT ASCENSION OF SUN FOR FOLLOWING DAY (0 HRS. UT ) 

C T ( 5 ) - DECLINATION OF SUN FOR LAUNCH DAY (0 HRS. UT ) 

C T ( 6 1 - DECLINATION OF SUN FOR FOLLOWING DAY (0 HRS. UT ) 

C T (7 )— UNIVERSAL TIME OF LAUNCH 
C TL- TIME AFTER LAUNCH 

C COMPUTES GREENWICH SIDERIAL TIME, RIGHT ASCENSION, HOUR ANGLE, 

C DECLINATION, OF THE SUN ANO X3,Y3,Z3 
C CONVERTS ALL IMPUTTED TIMES TO SECONDS 
CONST 1= (T(7)+TL)/(24. *3600 . ) 

TG= (T < 1)+C0NST1=M T ( 2 )— T( 1 ) )+T{ 7) +TL ) /240. 
ALPHAS=(T(3)+C0NST1*(T(4) — T (3)))/240. 

EPSILS= (ALPHAS-TG— XLONG) /CONV 
DELTAS=(T(5 ) + CON ST I- ( T 1 6 1— T ( 5 ) ) ) / ( CONV*240 . ) 

XL AT = XLAT/CONV 
X3=C0S(DELTAS)*SIN1EPSILS) 

Y3=— SIN <XLAT)*COS (EPSI LSI *COS( DELTAS >+COS(XLAT)*S IN (DELTAS) 
Z3=C0S ( XLAT >*COS ( EPS I LS ) *COS( DELTAS )+SIN(XLAT)*S IN (DELTAS) 
MAGNETOMETER SECTION 
CALCULATES X2,Y2,Z2 WITH AZMITH IMPUT 
CALL T I MCOR (TD,T1,PHIM1,N1,DPHIM) 

PHIM=DPHIM/CONV 
PIE=2.-PI 

IF ( PH IM-P IE 1409,409,410 
410 PHIM=PHIM— P I E 
409 D0TP1 = X3*=X4+Y3*Y4+Z3*Z4 

A1M=Z1**2+(X1*C0S(PHIM)+Y1*SIN(PHIM) )**2 
B1M=-2.*Z1*D0TP1 

C1M=DOTP1*=!=2-(X1*COS(PHIM)+Y1*SIN(PHIM) )**2 
A2M=A1M 

B2M=-2.*DOTPl=i'(Xl*COS(PHIM)+Yl*SIN(PHIM) ) 

C2H=D0TP l*-2— Z 1**2 
RAD21=B1M**2-4.*A1M*C1M 
IF (RAD21.GT .0 . ) GO TO 408 
NF0RM=30 
GO TO 998 

408 RAD 1= SORT (RAD21 ) 

SINE 1=(— B1M+RAD1 ) / ( 2.=* C A1H» 

S INE2= (— B1M— RADI )/(2.-AlM) 

RAD22=B2M**2-4.*A2M*C2M 
IF(RAD22.GT.O. ) GO TO 407 
NF0RM=31 
GO TO 998 

407 RAD2=SQRT (RAD22 ) 

COSE 1= (-B2M+RAD2 ) / ( 2. *A2M) 

C0SE2=(— B2M-RAD2)/(2.*A2M) 

IF (COSE 1 .GT.O. ) GO TO 302 
IF(COSE2.GT.O. ) GO TO 305 
NF0RM=22 
GO TO 998 

302 IF (C0SE2.LT .0. ) GO TO 304 
NF0RM=23 
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GO TO 998 

304 CGTHET=COSE 1 
GO TO 306 

305 COTHET=CO SE 2 

306 CH6CK=C0THET**2+SINE1**2 
CHECK1=C0THET**2+SINE2**2 
C0NE=1. -CHECK 

CONE 1= 1. -CHECK 1 
TEST = 0.1 
700 TEST = TEST*0.1 

IFtABS(CONE). LE.TEST) GO TO 710 
IF(ABS(C0NE1) .LE.TEST) GO TO 720 
NF0RM=2 1 
GO TO 998 

710 IF(ABS(C0NE1). LE.TEST) GO TO 700 

307 S I THET = S INE 1 
GO TO 310 

720 S ITHET = S I NE 2 

310 X2=C0THET*C0S(PHIM) 

Y2=C0THET*SIN<PHIM) 

Z 2 = S I THET 

A6=(X1*Y2-Y1*X2) /( 1. -DOT P 1**2) 

A7=(Yl*Z2-Zl*Y2)/< l.-D0TPl**2) 

B6= ( Z 1-DOTP 1*Z 2 ) / ( 1 • -DOTP 1**2 ) 

B7=(X 1-D0TP1*X2) /( l.-D0TPl**2) 

C1 = Y3*Z4-Z3*Y4 

C2=Z3*X4-X3*Z4 

C3=X3*Y4-Y3*X4 

D 1 = X3“D0TP 1*X4 

D2 = Y3-DOTP 1*Y 4 

D3=Z3-D0TP1*Z4 

SIZSIA( 1)=Z2*X4+A6*C1+B6*D1 

S I ZCO A( 1)=Z2*Y4+A6*C2+B6*D2 

COSZENI 1)=Z2*Z4+A6*C3+B6*D3 

C0SZEN(2)=X2*Z4+A7*C3+B7*D3 

IF(C0SZEN(2).LT.1.)G0 TO 599 

N FORM= 17 

GO TO 998 

599 SIZSI A(2)=X2*X4+A7*C1+B7*D1 
SIZC0A(2 )=X2*Y4+A7*C2+B7*D2 
DO 630 J= 1 , 2 

ZEN ( J ) =AT AN ( ( SQRT( l.-COSZEN( J)**2) ) /COS ZEN (J ) ) 

I F ( ZEN ( J) >600,601,601 

600 ZEN ( J ) =P I +ZEN ( J ) 

601 SINZEN{ J)=SIN(ZEN( J) ) 

SINAZM(J)=SIZSIA(J)/SINZEN( J) 

COSAZM( J)=SIZCOA( J)/SINZEN( J) 

TANAZM( J )=SINAZM( J)/COSAZM( J) 

AZM(J)=ATAN( TANAZM ( J ) ) 

IF(TANAZM( JJ.GT.O. ) GO TO 604 
I F ( SINAZM(J) .GT.O. ) GO TO 603 
AZM(J)=2.*PI+AZM(J) 

GO TO 610 

603 AZM ( J ) =P I +AZM ( J ) 

GO TO 610 

604 IF<$INAZM< J) .GT.O. ) GO TO 610 
AZM ( J ) =P I +AZ M ( J ) 

610 CODCOE( J)=COS(XLAT)*COSZEN< J)-SIN< XLAT ) *S INZEN ( J )*COSAZM< J) 
CODSIEt J)=SINZEN( J)*SINAZM( J) 

SINDEL( J)=COS(XLAT)*SINZEN< J)*COSAZM< J ) +S IN ( XL AT ) *COSZEN < J) 
DELTA ( J)=ATAN( ( SI NOEL ( J ) ) / ( SORT ( 1 . -S I NDEL ( J ) **2 ) ) ) 
COSDEL(J)=COS(DELTA(J) ) 

S INE P S ( J ) =C0 DSIE(J)/COSDEL(J) 

COSEPS(J)=CODCOE(J)/COSDEL(J) 

TANEPS(J)=SINEPS(J)/COSEPS(J) 

EPSIL(J)=ATAN(TANEPS(J) ) 


22 



IFCTANEPSC JI.GT.O.) GO TO 614 
IF(SINEPS(J) *GT .0 . ) GO TO 613 
EPSIL( J)=2.*PI+EPSIL( J) 

GO TO 620 

613 EPSILt J)=PI+EPSIL< J) 

GO TO 620 

614 IF(SINEPS(J)*GT*0.) GO TO 620 
EPSILC J)=PI+EPSIL< J) 

620 EPSILt J)=EPSIL( J)*CONV 
AZMf J)=AZM! J)*CONV 
ZENC J)=ZEN( J)*CONV 
OELT A ( J ) =DELTA( J ) *CONV 
ALPHA ( J)=TG+XLONG+£PSIL( J) 

IF ( ALPHA C J ) *LT* 360* ) GO TO 625 
ALPHA ( J )= ALPHA ( J )-360* 

625 XH{ J ) -ALPHA < J 5/15* 

IX { J)=XH(J) 

XXI=IX ( J ) 

XM ( J 5 = ( XH ( J } — XX I } *60 • 

NX ( J ) =XH ( J } 

XXN=NX ( J ) 

XS( J)=(XM{ J)-XXN)*60. 

630 CONTINUE 

IF (NPLOT *EQ*Q ) GO TO 631 
NOUNT =N0UNT+1 
CPLVAR < 1 , NOUNT) =TD 
CPLVAR { 2 » NOUNT )=AZM( 1 ) 

CPL VAR ( 3, NOUNT )=AZM( 2 ) 

CPLVAR (4, NOUNT l^ZENt 1 ) 

CPLVAR ( 5 , NOUNT } =ZEN ( 2 ) 

CPLVAR (6 t NOUNT 1=XH(1) 

CPLVAR (7, NOUNT >=XH( 2) 

CPLVAR ( 8 t NOUNT 5=DELTA C 1 5 
CPLVAR (9,N0UNT)=0ELTA( 2 ) 

631 WRITEC6,9) 

WRITE (6, 10) T IHE t ZEN ( 1 ) , AZH( 1 ) , ZEN( 2) , AZM< 2 ) 

WRITE (6,15) 

WRITE <6, 16) IX ( 1),NX< 1 ) , XS(15 , DELT A(l) , I X < 2 ) ,NX ( 2 ) , XS < 2 ) , DELTA { 2 ) 
WRITE (6,20) 

JUMP= JUMP+6 
GO TO 999 

998 WRITE ( 6 »NFORM ) TD 
WR I TE (6 , 20 ) 

JUMP= JUMP+3 

999 I F ( JUMP* LE *59 ) GO TO 98 
WRITEC6, 19) 

JUMP=0 
GO TO 98 

800 IF(NPLOT*EQ.O) GO TO 1000 
RE AD (5 * 3 ) ( AA( J ) » J= 1 , 12 ) 

READ ( 5 » 4 ) NPLTS 
CALL CPLTAP C 10 ) 

00 820 JPL0T=1,NPLTS 
RE AD ( 5 » 3 ) ( BB( J ) , J=l, 12 5 
RE AD I 5 * 3 ) (CC( J),J=1,12) 

READ (5,4) ( (NSCALEt J), J=l,5 ) ,NH,NV, JCPLT) 

RE AD ( 5 ,4 ) (IA,I0(N),N=1, JCPLT) 

READ (5,5) SBH ,SBV,X2,X1,Y2,Y1 
CALL CP LOTI (NSCALE,NH, $BH,NV, SBV) 

CALL CPLOT 2(X2,X1,Y2,Y1) 

KK = 1 

DO 810 JCPLT 1= 1 , JCPLT 
NO=IG< JCPLT1) 

CALL CPL0T3(BCDX(KK) , CPL VAR (I A, 1) , CPL VAR ( NO , 1 ) ,NOUNT ) 

810 KK=KK+1 

820 CALL CPL0T41 72, BB(1) ,72,CC(1) ,72, AA( 111 
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CALL ENDPLT 
1000 STOP 
END 

$ I BFTC TIMC LI ST, REF, DECK , M94 

SUBROUTINE TIMCOR ( TI M , Y , THET , MT , ANS ) 

C MT=NUMBER OF VALUES IN TABLE 

C ANS=VALUE OF THET CORRESPONDING TO TIM IN T VALUES 

DIMENSION Y ( 400 ), THET (400 ) 

AT IM=T IM-Y ( 1) 

IF (ATIM)3,4,5 

4 ANS=THET( 1) 

GO TO 50 

3 KT = 1 

99 FORMAT ( 38H VALUE IS OUTSIDE RANGE OF TABL E , T IME= , E 16 . 8 , 7H LIMIT=, 
XE 16 *8 ) 

97 WRITE ( 6 >99 ) TIM,Y(KT) 

GO TO 51 

5 BTIM=T IM-Y ( MT ) 

IF (BTIM) 6,7,8 

8 KT = MT 
GO TO 97 
7 ANS=THET ( MT ) 

GO TO 50 

6 DO 47 K=2,MT 
CT I M=T I M-Y ( K ) 

IF(CTIM) 45,46,47 

46 ANS=THET ( K ) 

GO TO 50 

45 DELT=Y ( K) -Y ( K — 1 ) 

TOEL=T IM-Y ( K-l ) 

R E L=TDEL/DE LT 

ANS=THET( K-l)+( THET ( K) -THET (K-l ) )*REL 

GO TO 50 

47 CONTINUE 
GO TO 50 

51 STOP 
50 RETURN 
END 

SIBFTC FIELD L I ST, R EF , DECK , M94 

SUBROUTINE F I ELD ( FLAT , FLONG , ALT , BT , BP , BR , B2 ) 

C DUMMY SUBROUTINE FOR SETTING UP MCILWAINS MAGNET 

DIMENSION V ( 3 , 3 ) 

V( 1,2) - ALT/6371.2 
V ( 2 , 2 ) = (90. -FLAT)/57. 2957795 
V ( 3 , 2 ) = FLONG/57. 2957795 
AER = V ( 1 , 2 ) 

SIT = ABS ( SIN (V( 2,2) ) > 

SSQ = SIT**2 

1002 I F (V ( 3,2) ) 1000, 1001, 1001 

1000 V ( 3 , 2 ) = V ( 3, 2 ) + 6.283185307 
GO TO 1002 

1001 CALL MAGNET (AER,SIT,V(3»2) ,BR,BT»BP,B2»V(2,2) ) 

CONTINUE 

RETURN 

END 

$ I B FTC MAGN L I ST , R EF , DECK , M94 

SUBROUTINE MAGNET ( R , S , PH I , BR , BTHET , BPH I , BB ,THET ) 

3 FORMAT (1H0,5X,F14.7,1X,F14.7,1X,F14.7,1X,F14.7) 

DIMENSION DP (49 ) ,P (49 ) , G ( 49 ) ,H(49) , CONST (49 ) ,A0R(7) ,SP(7),CP<7) 
IF(KIP)151, 150,150 
150 K I P = — 1 

DO 152 N=l,49 
G ( N ) =0 .0 
152 H(N)=O.0 

JENSEN AND CAIN COEFFICIENTS FOR 1960 (JUNE 1962) 

G ( I ) = G ( N , M ) AND H ( I ) = H ( N , M ) WHERE I = N + 7 : MM-1) 


MAGNTOOO 

MAGNT001 
MAGNT002 
MAGNT003 
MAGNT004 
MAGNTOO 5 
MAGN TOO 6 
MAGNTOO 7 
MAGNT008 
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G ( 2 ) = 3.04112050E-01 
G( 9 ) = 2* 147368 58 E~0 2 
G ( 3 ) = 2.40 35367 1E—0 2 
G( 10)=-5.12533379E-02 
G ( 17 )=— 1# 338 1 1969 E— 02 
G ( 4)=— 3.15178651E— 02 
G ( 11)= 6. 21300906E— 02 
G ( 18 )= -2* 4898 1333 E— 02 
G ( 25 ) =— 6 • 49565905 E— 03 
G ( 5 3=— 4. 17943639E— 02 
G ( 12 ) =— 4 .5298 3660 E— 02 
G ( 19 >=-2. 17947447 E— 02 
G ( 26 ) = 7. 008 2540 5 E— 03 
G ( 33 ) =— 2 .0439 5562 E— 03 
G ( 6 )= 1.62556271E-02 
G { 13 )-— 3.44067606E-02 
G ( 20 )=-1.94470026E~02 
G(27)=— 6.08211374E-04 
G ( 34 ) = 2 . 775 33 549 E— 03 
G ( 4 1 ) = 6 .96802467 E— 04 
G ( 7 )=— 1.952317 36 E— 02 
G ( 14 )=— 4.85326147E— 03 
G ( 21 ) = 3.21172 428E— 03 
G { 28 ) = 2. 14128828E-02 
G f 35 ) = 1 .0505 1275 E— 03 
G {42 ) = 2 • 26829448 E— 04 
G ( 49 ) = 1 • 1147 1358 E— 03 
H ( 9 ) = — 5 .7989050 IE— 02 
H ( 10 ) = 3. 312407 14E— 02 
H( 17) =-1.5789 3822E-03 
H ( 1 1 ) = 1. 48696943E— 02 
H ( 18 ) =— 4 .07490158E— 0 3 
H C 25 ) = 2. 10318235E— 04 
H ( 1 2 ) =- 1 . 18 245456E— 02 
H ( 19 ) = 1 .00057732 E— 02 
H { 26 ) = 4 • 3038 086 3 E— 04 
H ( 33 )= 1 • 38 50 3490 E —03 
H( 13)=— 7.95897466E-04 
H ( 20 ) =— 2 • 000440 21 E— 03 
H { 27 ) = 4.59718859 E— 0 3 
H { 34 ) = 2 • 4206 3078 E—0 3 
H ( 4 1 ) =-l • 2 1806522 E—0 3 
H( 14) =-5.758 30293E-03 
H< 21 )=-8. 7346 1401 E-03 
H ( 28 ) =— 3. 40604073E— 0 3 
H(35 )=-l. 18 16 2456E— 04 
H(42)=-l. 1 16230 13E— 03 
H ( 49) =-3.2483 189 IE-04 
P(l)=1.0 
DP { 1 ) =0 • 0 
SP( 1)=0.0 
CP( 1)=1.0 
CONST ( 9 ) =0 *0 
CONST (16) =0.0 
DO 80 N = 3f 7 
FN=N 

DO 80 M= 1 , N 
FM = M 

I=N+7*CM-1) 

80 CONST (I ) = ( (FN-2.0 )**2-<FM-1.0)**2)/( ( FN+FN-3 .0 ) * ( FN+FN-5 . 0 ) ) 

151 C=SQRT (ABS (1.0-$*S)) 

IF(THET-1. 570796327) 154,154,156 

156 C=— C 

154 AR=1./(1.+R) 

155 SP ( 2 ) =SIN (PHI) 

CP ( 2 ) =COS (PHI) 


MAGNT009 
MAGNT010 
HAGNT011 
MAGNT012 
NAGNT013 
MAGNT014 
MAGNT015 
MAGNT016 
MAGNT017 
MAGNTO 18 
MAGNTO 19 
MAGNT020 
MAGNT021 
HAGNT022 
MAGNT023 
MAGNT024 
MAGNT025 
MAGNT026 
MAGNT027 
MAGNT028 
MAGNT029 
MAGNT030 
MAGNT031 
MAGNT032 
MAGNT033 
HAGNT034 
MAGNT035 
HAGNT036 
MAGNT037 
MAGNT038 
MAGNT039 
MAGNT040 
MAGNT041 
MAGNT042 
MAGNT043 
MAGNT044 
MAGNT045 
MAGNT046 
MAGNT047 
MAGNT048 
HAGNT049 
MAGNT050 
MAGNT051 
MAGNT052 
MAGNT053 
MAGNT054 
MAGNT055 
MAGNT056 
MAGNT057 
MAGNT058 
MAGNT059 
MAGNT060 
MAGNT061 
MAGNT062 
MAGNT063 
MAGNT064 
MAGNT065 
MAGNT066 
MAGNT067 
MAGNT068 
MAGNT069 
MAGNT070 
MAGNT071 
MAGNT072 
MAGNT073 
MAGNT074 
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A0R< 1)=AR*AR 

‘ MAGNT075 


AOR ( 2 ) = AR*AOR ( 1 ) 

MAGNT076 


DO 90 M = 3* 7 

MAGNT077 


N=M- 1 

MAGNT078 


SP(M)=SP(2)*CP(N) +CP ( 2 ) *SP ( N ) 

MAGNT079 


CP (M)=CP(2)*CP(N)-SP(2)*SP(N) 

MAGNT080 

90 

AOR ( M ) = AR*AOR ( N ) 

MAGNT08 1 


BR = 0 *0 

MAGNT08 2 


BTHET=0 *0 

MAGNT08 3 


BPH I =0 *0 

MAGNT084 


DO 32 N = 2, 7 

MAGNT08 5 


FN=N 

MAGNT08 6 


SUMR=0.0 

MAGNT087 


SUMT=0.0 

MAGNT088 


SUMP=0 • 0 

MAGNT089 


DO 33 M=1,N 

MAGNT090 


I F ( N — M )87, 88,87 

MAGNT091 

88 

I =8*N-7 

MAGNT092 


L = I —8 

MAGNT093 


P(I )=S*P(L> 

MAGNT094 


DP € I )=S*DP(L)+C*P(L) 

MAGNT095 


GO TO 89 

MAGNT096 

87 

I=N+7*(M-1> 

MAGNT097 


J = I -1 

MAGNT098 


K = I-2 

MAGNT099 


P ( I )=C*P( J) -CONST (I ) *P ( K ) 

MAGNT100 


DP (I )=C*DP( J)-S*P( J) -CONST ( I )*DP(K) 

MAGNT101 

89 

FM=M— 1 

MAGNT102 


TS=G ( I )*CP(M)+H(I )*SP(M) 

MAGNT103 


SUMR=SUMR+P ( I ) *T S 

MAGNT104 


SUMT=SUMT+DP ( I )*TS 

MAGNT105 

33 

SUMP=SUMP+FM*P( I ) ^ ( -G ( I ) *SP ( M ) +H ( I ) *CP ( M ) ) 

MAGNT106 


BR=BR+AOR { N ) *FN*SUMR 

MAGNT107 


BTHET = BTHET-AOR ( N ) *SUMT 

MAGNT108 

32 

BPHI =BPHI -AOR { N ) *SUMP 

MAGNT109 


BPH I = BPH I / S 

MAGNT110 


BB = SQRT ( BR**2+BTHET**2+BPHI**2) 
WRITE (6,3) BR, BTHET, BPHI ,BB 

MAGNT111 


RETURN 

MAGNT112 


END 

MAGNT113 
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8.0 RESULTS AND CONCLUSIONS 


The orientation of a rocket can be obtained in horizon or space-fixed coordinates by obtaining 
the spatial location of the solar and magnetic field vectors with respect to both a rocket- axis and 
horizon coordinate system, with the use of only a lateral magnetometer and solar sensors. 

The accuracy of the procedure is increased as the spin rate of the vehicle increases or becomes 
closer to a constant value. It is also difficult to obtain a solution if the solar and magnetic field 
vector are nearly coincident in space at a given time since two distinct vector are required in the 
procedure. 

Another limitation of the method is the accuracy of the various sources of data which are used. 

(1) The location of the magnetic vector in horizon coordinates is obtained from a 48 term expansion 
in spherical harmonics. This expansion becomes slightly less accurate at very high altitudes. Care 
must be taken to insure that the most recent magnetic data has been used in the expansion. (2) The 
location of the solar vector in horizon coordinates obtained from reference (2) is quite accurately 
known and any resulting error is negligible; however, both (1) and (2) require accurate radar data. 
The accuracy of this data depends of course on the particular system and smoothing procedure used. 
(3) The accuracy of the solar sensors which determines the location of the solar vector with re- 
spect to the rocket-axis coordinate system is dependent on the particular type of sensor used; how- 
ever, the majority of the sensors use the sun as a point source which results in a minimum error 
of approximately one-half of a degree which corresponds to the width of the solar disc. This error 
can be reduced somewhat by smoothing the sensor data timewise. (4) The accuracy of the mag- 
netometer which determines the orientation of the magnetic field vector in the rocket-axis system 
is a function of the precision of the instrument used. The accuracy of the method is not affected 
by any magnetic source in the rocket payload since this effects only the absolute amplitude of the 
lateral magnetometer output. The method requires only a knowledge of the timewise location of the 
maxima and minima of the output. 

It should also be noted that the method assumes that the sun is within view of at least one 
sensor. It is usually difficult to have a sensor mounted at the tail of the rocket; therefore, no 
can be obtained if the tail is facing the sun or the flight occurs during a time when the sun is 
visible. 

9.0 REFERENCES 

(1) National Almanac Office, U.S. Naval Observatory: The American Ephemeris and National 
Almanac for the Year 1963. United States Government Printing Office, 1961. 

(2) Cain, J. C., Hendrick, S., Daniels, W. E., Jensen, D. C.: Computation of the Main Geomagnetic 
Field from Spherical Harmonic Expansions, NASA X-61 1-64-316, 1964. 

(3) Phillips, H. B.: Vector Analysis. John Wiley & Sons, Inc., cl933. 

(4) Thomson, W. T.: Introduction to Space Dynamics. John Wiley & Sons, Inc. 1961. 


solar 

data 

not 


27 


APPENDIX A 
Illustrations 


ROCKF.T LONGITUDINAL AXIS 



/ 

X 

(SENSOR AXIS) 


Figure 1 -Location of Solar Vector in Rocket-Axis System using Case I Solar Sensor (Lateral Sensor) 


ROCKET LONGITUDINAL AXIS (SENSOR AXIS) 



Figure 2-Locotion of Solor Vector in Rocket-Axis System using Case II Solar Sensor (Longitudinal Sensor). 



NORTH CELESTIAL POLE 



ROCKET POSITION 


Figure 4-Location of o Vector in the Equatorial Plane. 
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VERTICAL AXIS 



NORTH AXIS 


Figure 6-Location of the Magnetic Field Vector in the Horizon Coordinate System. 
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iTIAL POLE 



VERNAL EQUINOX 

Figure 8-Location of a Vector in Space-Fixed Coordinate System. 
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APPENDIX B 


Sample Results 

The procedure was applied to an Aerobee rocket and representative time histories of Azimuth 
and Zenith angles and Right Ascension and Declination for both the Experiment and Rocket Longi- 
tudinal axes are shown. 

Aerobee 4.13 GP-GT was launched September 26, 1963. The rocket was de-spun from approxi- 
mately 2.5 rps after burnout at which time an attitude control system was to have controlled the 
vehicle orientation. This system, however, did not function properly and the vehicle experienced a 
reverse roll rate of about 0.05 rps. A coning motion resulted with the coning and roll periods being 
approximately equal. The results presented describe the orientation during the reverse roll (170- 
300 seconds). 
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